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Lecture  8 

COUNTERFLOW  DIFFUSION  FLAMES 

The  fundamental  characteristic  of  diffusion,  flames  is  that  the  two 
reactants,  fuel  and  oxidizer,  are  supplied  in  different  parts  of  the 
combustion  field,  so  that  they  must  come  together  and  mix  by  diffusion 
before  reaction  can  take  place.  Count erf loving  streams  provide  one 
method  of  bringing  then  together;  the  resulting  diffusion  flane ,  whose  * 

s» 

main  properties  were  established  by  Lina n,  is  the  subject  of  this  lecture. 

1.  Basic  Equations 

Consider  the  combustion  field  sketched  in  figure  1.  A  stream  of  gas 
containing  the  oxidant  Y1  =  X  flows  to  the  left  and  iapinges  on  a  stream 
containing  the  fuel  Y 2  =  Y  that  flows  to  the  right,  forming  a  stagnation 
point  at  the  origin.  The  flow  field  is 

(u,v)  *  2(-x,y)  (1) 

•under  the  constant -density  approximation ;  here  a  constant  of  proportionality 

s/2,  where  e  is  the  straining  rate,  has  been  absorbed  into  the  length 

'unit  (which  is  used  to  define  M  ).  For  such  a  flow,  it  is  possible  for  the 

r 

combustion  field  to  be  stratified  in  the  y-direction,  with  a  flat  flame 
sheet  at  x  =  x#.  The  temperature  and  mass  fractions  are  then  functions 
of  x  and  t  alone  satisfying,  for  unit  Lewis  numbers, 

—8  /T 

L(T)  *  -2L(X)  =  -21(f)  *  with  fl  »  VXYe  ;  (2) 

1  =  3/3t-2x  3/3x  -  32/3x2  (3) 

and  the  boundary  conditions  are 

-8.1- 


—8 . 2— 


2  *  T,_,  X  0,  Y  ->  Yf  as  x  ->  X  -  Xf ,  Y  -  0  as  X-+-.  (U) 

The  coefficient  17  is  proportional  to  1/e,  so  that  an  increase  in  the 
straining  rate  causes  a  decrease  in  P.  It  is  the  response  of  the 
combustion  field  to  variations  in  P  that  is  of  principal  interest. 

Since  "both  Levis  numbers  have  been  taken  equal  to  1,  there  are  two 
Shvab-Zeldovich  variables  (section  2.2),  namely 

T+2X  *  T++Z_,  T+2Y  =  T++Z+,  (5) 

where 


?+  =  J[Tf +erfc(-x)+Tf_erfc(x) ] ,  Z_  =  Xferfc(-x),  Z+  *  Yferfc(x).  (6) 


3oth  these  variables  are  annihilated  by  1  and  satisfy  the  boundary 
conditions  (U ) .  We  are  left  with  a  problem  for  the  temperature  alone, 
when  X'  and  Y  are  hereby  suppressed  in  favor  of  T.  Note  that  no 
assumption  of  steadiness  has  been  made,  but  these  relations  can  only  be 
applied  to  unsteady  problems  for  which  the  initial  conditions  satisfy 
the  relations.  We  shall,  however,  be  concerned  solely  with  steady  solutions 
until  section  5- 
When 

_q  /m 

(7) 


d.-6/t  -  o 


there  is  no  chemical  reaction,  even  if  reactants  are  present,  and  the 
combustion  is  said  to  be  frozen.  There  are  two  ways  in  which  this  can 
be  brought  about: 

(i)  P  0,  the  small  Daakohler -number  limit; 

(ii)  P  e®^1*  as  0  -*■  «•  with  T  <  T„,  a  creature  of  activation- 

•  energy  asymptotics. 

In  case  (i)  the  entire  combustion  field  is  frozen,  whereas  in  case  (ii) 
only  that  portion  where  T  falls  below  T#  is  frozen. 


-8.3- 


Wher. 


(8) 


there  is  nothing  to  balance  an  infinite  reaction  rate  outside  vanishingly 
thin  layers  (spatial  or  temporal),  and  so 

XY  -►  0.  (9) 

4 

At  least  one  of  the  reactants  is  not  present,  i.e.  there  is  equilibrium. 
Here,  as  for  frozen  combustion,  there  is  no  reaction,  but  for  quite  a 
different  reason. 

There  are  two  ways  in  vhich  equilibrium  can  be  achieved: 

(i)  V  -*■  ■*•,  the  large  Damkbhler  number  limit; 

0/T„ 

(ii)  V  -  e  as  8  -►  ®  with  T  >  T#,  another  creature  of  activation- 

energy  asymptotics. 

In  case  (i)  the  entire  combustion  field  is  in  equilibrium,  but  in  case  (ii) 
only  that  portion  where  T  rises  above  T#  is  in  equilibrium-  Because 
of  the  nature  of  this  limit,  there  is  always  the  possibility  of  a  reaction 
zone,  known  as  a  Burke-Schuaann  flame  sheet,  existing  In  the  middle  of  an 
equilibrium  region;  that  cannot  happen  In  a  frozen  region. 

Simple  analytical  treatments  are  possible  In  these  two  limits, 
frozen  and  equilibrium.  So  elementary  Is  the  first  that  we  shall  forego 
discussion  of  it,  concentrating  instead  on  the  second,  which  is  called 
the  3urke-Schumanr.  (equilibrium)  limit  after  an  early,  basic  combustion 
problem  they  considered  (see  section  10.2).  Large  values  of  V  are 
easily  obtained  in  practice,  the  reason  why  the  limit  was  originally 
introduced  by  Burke  and  Schumann. 

To  the  left  of  the  flame  sheet  at  x  =  x„  there  is  no  oxidant,  and 
to  the  right  no  fuel.  The  Shvab-Zeldovich  relations  (5)  therefore  give 
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HVZJ 

for  x  >  x*,  (10) 

J 

and  continuity  across  the  flame  sheet  then  requires  Z_  =  Z+  there,  i.e. 

erf(x*)  =  (VXf'/(-f+Xf)  (H) 

which  determines  x„.  The  flaae  sheet  lies  in  the  leaner  stream:  x„  >  0 

accordingly  as  Xf  *  Yf.  The  flame  temperature  is 

T#  *  {Tf_xf +  Vf  +  2XfYf)/(xf+yf)i  (12) 

its  structure  (35),  (3o)  will  be  derived  in  section  3  (for  9  —). 

2.  The  S-Shared  Burning  Response 

The  usual  way  of  characterizing  the  solution  is  to  plot  variations  in 
some  significant  parameter,  such  as  the  maximum  temperature,  with  the 
Da-kShler  number.  If  Tf+ ,  Xf,  Yf  are  such  that  the  combustion  generates 
a  heat  flux  to  both  far  fields,  this  response  is  S-shaped  in  the  limit 
8  ■*  -(figure  2).  Certain  physical  conclusions  can  then  be  drawn. 

If  the  system  is  in  a  state  corresponding  to  a  point  on  the  lower 
branch,  and  V  is  slowly  increased,  the  solution  can  be  expected  to 
change  smoothly  until  the  point  I  is  reached.  Rapid  transition  to  the 
upper  branch  will  then  presumably  occur,  corresponding  to  ignition.  A 
subsequent  3low  decrease  in  V  is  likewise  anticipated  to  produce  a 
smooth  decrease  in  burning  rate  until  extinction  occurs  at  E. 

If  one  of  the  far  fields  loses  heat,  the  response  is  monotonic,  so  that 
the  phenomena  of  ignition  and  extinction  are  absent.  For  high  activation 
energy  the  transition  from  a  monotonic  to  an  S-shaped  response  occurs 
when  the  tenrcerature  gradient  on  one  side  of  the  flsme  sheet  is  small,  a 
case  that  is  not  difficult  to  analyze. 


' 

z 

0 

"  ,  X  = 

K 

II 

Hz  -z  ) 
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Assume,  without  loss  of  generality ,  that 


V  ”  Tf-; 


then  the  sr.all  temperature  gradient  rust  be  on  the  right  of  the  flame  sheet, 
i.e.  hotter  side.  If  it  were  on  the  colder  side,  the  flame  temperature 
(12)  would  be  close  to  ,  which  implies  (in  the  limit,  of  zero  gradient) 


Xf  *  J(Tf_-Tf+)  <  0,  (1U) 

ar>  impossibility.  To  ensure  that  Tf+  is  dose  to  the  flame  temperature, 
ar.d  hence  that  the  temperature  gradient  is  small  on  the  hotter  side,  set 


2 Yf  =  Tf+-Tf_  +  k/9  with  k  »  const.  (15) 

Ir.  seeking  an  asymptotic  solution  as  e  ■+■  «*  we  shall  assume  that  equilibrium 
prerails  for  x  >  x#  even  though  T  does  not  rise  above  T#  by  an  0(1) 
amount  there ,  and  check  a  posteriori  that  the  solution  thereby  constructed 
is  self-consistent. 

In  *.’iew  of  the  assumption  (15),  the  Shvab-Zeldovich  relations  (5b) 
becomes 

T  +  2Y  -  Tf+  +(k/29)erfc(x)  (l6) 


r>  a 

—  *>  *—**'''■  9 


T  *  7  +  (k/29 )erfc(x )  for  x  >  x#,  (17) 

since  Y  =  0  there.  To  complete  the  description  of  the  combustion  field 
one side  the  reaction  zone,  we  need  the  temperature  in  the  frozen  region 
ahead  of  the  flame  sheet,  i.e.  the  linear  combination  of  1  and  erf(x)  that 
takes  on  the  values  T^,  at  x  =  -»  and  T#  at  x  *  x#>  clearly 
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T  =  (T,  [erfc(-x#)-erfc(-x)l  +  T#erfc(-x)}/erfc(-x#)  i'or  x  <  x#  .  (l8) 

This,  like  the  result  (17),  is  correct  to  a ry  order  in  3  provided  T# 
is  determined  to  the  same  order.  We  shall  only  need  leading-order 
accuracy  in  the  result  (18),  so  that  taking 

T„  *  Tf+  (19) 

is  good  enough.  Determination  of  the  non-exp amiable  x#  comes  from  analysis 
of  the  reaction  zone,  for  which  the  leading-order  result 

X#  »  c(Tf_-T_+)erfc(x#)+jXferfc(-x#) ,  (20) 

a  consequence  of  the  Shvab-Zeldovich  relation  (5a),  is  needed. 

The  appropriate  variable  in  the  reaction  zone  is 

£  =  8 ( x-x# ) ,  (21) 

so  that  coefficients  in  the  layer  expansion 

T  »  T  -  e"1^  +...  with  <J>  =  (l/TL  (22) 

f+  f+  1 

are  considered  to  be  functions  of  £.  The  Shvab-Zeldovich  relation  (5b) 
gives 

Y  =  jS'‘1T2(d-<J’ *)+•••  with  =  -k  erfc(x„)/2T2,  (23) 

f+  *  f+ 

so  that  the  structure  equation  is 

oo'  a  -O'?  o 

d"d/di  =  £Kd-$#)e  9  with  V  =  PX„e  *720  .  (2fc) 

rrote  that  V  =  0(1)  implies  Ve  '  *  0(9  )  for  x  >  x„:  there  must  be 
the  equilibrium  Y  *  0  on  the  right  of  the  flame  sheet  since  otherwise 
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the  reaction  term  ir.  the  Y-equation  could  r.ot  be  balanced  there.  In  other 
vords,  our  solution  is  self-cansister.t . 

Equation  (2h)  is  precisely  the  structure  equation  for  the  prenixed 
flames  discussed  ir.  section  2.h.  Since  the  gradient  dT/dx  vanishes 
or.  the  right  of  the  flame  sheet ,  it  determines  the  gradient  on  the  left 
(as  was  explained  in  section  2.5)*  But  the  latter  is  known  in  terms 

« 

of  xH  from  the  expression  (18),  so  the  result  is  an  equation  for  x#, 
namely 


2(W 

7tt 


-X, 


erfcT-xTT 


2  «-•  -O*12 
~  0  /  2? 


6e 


f+ 


•) 


When  the  definition  (20)  is  used  to  eliminate  X», 

equations  (23)  and  (25)  give  (representing,  when  it  is  negative, 

the  maximum  temperature  in  the  combustion  field)  and  V  as  functions  of 
x#,  i.e.  the  required  relation  between  the  maximum  temperature  and  the 
DamkShler  number. 

The  corresponding  response  curve  is  shown  in  figure  3  for  several 
positive  values  of  the  constant  k;  when  k  is  non-positive,  the 
maximum  temperature  is  Tf+  and  not  T^.+-0  ^T2+b*.  For  k  sufficiently 
small,  the  response  is  monotonic;  otherwise  it  is  S-shaped.  Responses 
in  the  shape  of  an  S  appear  to  be  associated  with  a  flux  of  heat 
away  from  the  flame  sheet  on  both  sides  ,  but  this  has  never  been  proved. 


3.  General  Extinction  Analysis. 

We  turn  now  to  the  question  of  extinction  in  general,  i.e.  when  there 
is  an  0(1)  heat  flux  away  from  the  flame  sheet  in  both  directions.  It  is 
found  that  the  flame  temperature  on  the  whole  upper  branch  of  the  S-response 
differs  by  0(8-^)  from  the  Burke -Schumann  value  (12)  at  its  end;  our 
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extinction  analysis  vill  use  this  fact.  3ecause  of  the  0(1)  droo  ia 
temperature  avay  from  the  flame  sheet,  the  combustion.'  field  or.  each  side 
is  now  frozen-,  nevertheless,  the  leadir.g-order  solution  outside  the 
reaction  zone  is  identical  to  the  equilibrium  solution  constructed  in.  the 
limit  V  -*■  ®  (section  1).  Thus,  the  results  (10 ) ,  vhich  follow  from 
equilibrium  and  the  Shvab-Zelaovich  relations,  hold  for  ,XQ , -q-  Finally, 
before  considering  the  flame-sheet  structure,  we  define  x#  and  T#  by* 
the  formulas  (ll)  and  (12). 

In  the  reaction  zone,  the  variable  (2l)  and  expansion  (22)  still  applv, 
while  the  Shvab-Zeldovich  relations  gives 

X  =  50_1tJ(<|>+A5  )+...,  Y  =  J8"1rj(*+  35)+...,  (26) 

where 

2  2 

A  =  (Tf+-Tf_+2Xf)e  *MrT#>0,  B  =  (Tf+-Tf_-2Yf  )e  */At2  <  0;  (27) 


the  signs  of  A  and  3  follow  from  the  requirement  that  the  flame 
temperature  be  larger  than  the  temperatures  at  both  -»  and  +®.  The 
temperature  equation  therefore  reduces  to 


d2<f>/d52  =  $(4>+A5)($+B5)e-<*  with  V  =  T2  V, 


-9/T 


Ve3, 


(28) 


under  the  boundary  conditions 


d<>/d£  =  -A+...  as  5  dd/d£  =  -B+...  as  5  ■+  +°°,  (29) 

vhich  come  from  matching  the  leading  terms  in  the  expansions  (26)  with 
Xq  *  0  for  x  <  x*  and  *  0  for  x  >  . 


(30) 
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The  problem  (28),  (29)  determines  the  minimum  value  (i.e.  the  maxiumum 

memoeratu re  in  the  combustion  field)  as  a  function  of  0.  It  can  be  reduced  to 
canonical  form  by  writing 

5  =  (J-'-yS*  l  =  (A-B)c/2  with  y  =  (.-+3)/(A-B);  (31 ) 

vs  find 

d26/di2  =  5  (^2-52)e"^^  vith  O'  =  UP/(A-B)2  (32) 

e  “ 

and 

db/dC  =  -1+...  as  £*►-»,  d$/d£  =  1  +...  as  £  -*■  +«.  (33) 

Various  responses  determined  by  numerical  integration  of  the  problem  (32), 

(33)  are  shown  in  figure  U;  these  C-shaped  curves  correspond  to  the  neigh¬ 
borhood  of  the  point  E  in  figure  2.  Note  that  A  >  0,  B  <  0  imply 

-1  <  Y  <  1.  (34) 

Certain  limiting  forms  of  this  problem  as  t>e  -*■“  are  of  interest. 

(i)  The  Burke-Schumann  limit 

d2$/dC2  =  U2-C2),  d$/d£  =  +1+...  as  I  -*■  +”  (35) 

is  obtained  by  setting 

5=0  "l/30  ,  I  =  (36) 

e  e 

This  is  the  structure  problem  for  the  equilibrium  solution  discussed  at  the 
end  of  section  1;  it  applies  far  to  the  right  on  the  upper  branch  in 
figure  U .  Existence  and  uniqueness  have  recently  been  proved  by  Holmes. 

(ii)  The  so-called  preraixed- flame  limit  results  on  defining  a  small 
parameter  e  by 
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(9=0).  They  are  of  interest  because  every  point  on  the  middle  branch 
of  the  S  in  figure  2  that  is  not  too  close  to  E  or  1  corresponds 
to  a  solution  for  which  the  f lane-sheet  structure  has  one  of  these 
three  forms.  As  an  example  ve  shall  now  demonstrate  the  general 
applicability  of  case  (iii). 

4.  Partial-Burning  Branch  , 

This  part  of  the  S  is  characterized  by  0(1)  values  of  X*  and  Y„, 
hence  the  tern  partial-burning.  It  is  convenient  to  prescribe  the  flame 
temperature  (also  the  maximum  temperature)  and  calculate  V,  rather 
than  vice-versa. 

On  either  side  of  the  flame  sheet  the  leading-order  temperature  is 

Tq  *  Tf_  +  (T„-Tf_)erfc(+x)/erfc(+x1j)  for  x  $  x#.  (44) 

In  order  to  determine  x#,  we  anticipate  a  conclusion  to  be  drawn  from 
the  flame-sheet  structure,  to  wit 

+  §  (*.*»  »  0;  (U5) 

this  gives 

erfc(x*)/erfc(-x*)  =  (T#-Tf+)  (T#-Tf_) .  (46) 

The  Shvab-Zeldovich  relations  (5)  then  determine  X*  and  Y#. 

In  the  reaction  zone,  the  variable  (21)  and  expansion  (22)  still  apply, 
so  that  the  structure  equation  is 

-0  /T 

d2$/dC2  =  Ve~ *  with  V  =  PX#Y#e  */T^8,  (47) 

which  is  equivalent  to  the  equation  (43).  A  single  integration  shows  that 
the  derivative  d$/dx  takes  equal  but  opposite  values  at  £  *  ±»,  the  origin 
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of  the  relation  (1*5).  Mo  information  is  obtained  about  P,  ho -ever ,  except 

a/r* 

that  it  is  proportional  to  e  and  therefore  a  rapidly  decreasing  function 
of  T#;  to  determine  it  requires  an  examination  of  higher-order  terns. 

This  structure  is  appropriate  for  the  lover  part  of  the  niddle  branch. 

As  T#  increases  on  moving  up  the  branch,  one  of  the  mass  fractions  X# 
or  Y#  decreases  to  zero;  therefore,  a  different  structure  takes  over, 

_i  *  - 

characterized  by  0(9  *)  values  of  one  of  the  mass  fractions  but  still  0(1) 

values  of  the  other.  One  or  other  of  the  equations  (39)  then  governs; 

8/T, 

again  V  is  proportional  to  e  .  Further  increase  in  T#  towards  the 
Burke-Schunann  value  (12)  causes  the  remaining  0(1)  mass  fraction  to 
decrease  until  the  structure  (32),  (33)  is  attained.  We  forego  further 
discussion  of  the  middle  branch  since  it  corresponds  to  unstable  solutions, 
our  final  topic. 

5.  Stability. 

The  middle  branch  has  long  been  believed  to  be  unstable,  but  only 
recently  has  the  natter  been  confirmed  mathematically.  To  do  so,  the 
problem  must  be  examined  on  a  time  scale  that  is  relevant  to  the  reaction 
zone,  i.e.  using  a  fast  time 

t  =  e2t.  (U8) 

The  governing  equations  (2)  show  that 

3T/3t  =  9X/3r  =  3Y/3t  *  0  (1*9) 

everywhere  on  the  two  sides  of  the  flame  sheet:  T,  X,  and  Y  are  described 
to  all  orders  by  the  steady  state . 

On  the  other  hand,  in  the  reaction  zone  the  time  derivatives  are  as 
inoortant  as  the  diffusion  terms,  so  that  equation  (28)  becomes 
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325/3S2-3<!>/3t  =  P(d+A5)($+  (50) 

Correspondingly,  equation  (32)  is  replaced  by 

32$/352-3$/3^  ■  V^l2-l2U~^a  with  t  =  (A-3)2tA.  (51) 

Infinitesimal  disturbances  e**^  of  the  steady  state  vhich  is  the 

solution  of  the  problem  (32),  are  then  described  by  the  eigenvalue 
problem 

+  [V(5)-X]^  =  0,  ^  -  0  as  I  +  ±«,  _  (52) 

where 

-  -2  “Vy5 

V  *  t?e(«0-2«0-5  }e  °  .  (53) 

(Attention  is  restricted  to  disturbances  satisfying  the  Shvab-Zeldovich 
relations. ) 

This  problem  has  been  treated  by  Taliaferro,  Buckmaster  &  Nachman 
(1983)  vho  show,  in  particular,  that  the  transition  from  stable  solutions 
on  the  upper  branch  to  unstable  solutions  on  the  middle  branch  occurs 
exactly  at  the  static  extinction  point  E  in  figure  2.  (The  rigorous 
part  of  the  analysis  is  due  to  Taliaferro.)  Stability  on  the  upper  branch 
is  typified  by  the  Burke-Schumann  limit  (35),  for  which  the  corresponding 
eigenvalue  problem  is 

-  (2^+1)^  =  0  with  A  =  P~  JX,  <t>1  -*•  0  as  £  '►  (5*0 

Instability  on  the  middle  branch  is  exemplified  by  the  premixed-flame  and 
partial-burning  limits  (39)  and  A3).  For  the  former  the  eigenvalue 


oroblem  is 
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<>1  -  [2(li£-60)e  0  +X]<»1  *  0 

under  the  boundary  conditions  (5^b);  for  the  latter  it  is 

+  <e  -  0 


(55) 


(56) 


under  the  sane  boundary  conditions. 

Standard  treatments  shov  that  the  spectrum  of  X  is  negative  if 
is  everywhere  positive,  and  that  is  ensured  by  the  property  >  ||j  of 
the  steady  state.  Peters  was  the  first  to  show,  albeit  nunerically,  that 
the  problems  (55)  have  eigenvalues  with  positive  real  part.  (Peters 
writes  equations  in  a  slightly  different  way  from  ours;  his  m  is  equivalent 
to  2(l-|y|)  and  so  is  positive.)  Matalon  and  Ludford  noted,  in  a  different 
context,  that  the  problem  (56)  has  a  positive  eigenvalue.  Indeed,  the 
steady  state  is 

»  £n(2  cosh2J(£-Cm)],  (57) 

where  |  is  an  integration  constant  giving  the  location  of  the  mav-timim 
in 

temperature,  so  that 

X  *  l,  *  sech  lU-lj  (58) 

are  seen  to  be  an  eigenvalue  and  its  corresponding  eigenfunction. 

6.  The  Ignition  Point. 

The  neighborhood  of  the  point  I  in  figure  2  can  be  analyzed  by 
considering  perturbations  of  the  frozen  solution  that  raise  the  maximum 
temperature  by  an  0(0”^)  amount.  A  similar  analysis  in  the  context  of 
spherical  diffusion  flames  will  be  given  in  the  next  lecture,  so  we  omit 
the  discussion  for  counterflow  flames.  The  stability  of  the  lower  branch 
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has  not  yet  bean  examined,  but  the  conclusions  rust  agree  vith  the 
instability  result  for  the  partial -burning  portion  of  the  middle  branch. 
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Captions 


8.1  Notation  for  the  counterflow  diffusion  flame. 

8.2  Ignition  and  extinction  for  S-shaped  response. 

8.3  Steady-state  responses  when  the  temperature  gradient 
in  positive  for  x  >  0  and  small  for  x  <  0  ,  the 
latter  being  representative  by  K.  Drawn  for 

8.4  Extinction  curves.  For  |y|  ^  1  no  turning  point  is 
found  numerically,  a  result  that  still  lacks  formal 
proof. 
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